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Abstract 



One of the largest H a rotation curve data bases of 
spiral galaxies currently available is that provided 



by Persic & Salucci, Qastro-ph/9502091| ) hereafter PS 
1995, which has been derived by them from unre- 
duced rotation curve data of 965 southern sky spirals 
obtained by Mathewson, Ford & Buchhorn, hereafter 
MFB 1992. Of the original sample of 965 galaxies, the 
observations on 900 were considered by PS 1995 to 
be good enough for rotation curve studies, and the 
present analysis concerns itself with these 900 rota- 
tion curves. 

The analysis is performed within the context of the 
basic hypothesis that the phenomenology of rotation 
curves in the optical disc (that is, away from the dy- 
namical effects of the bulge) can be systematically 
described in terms of a general power-law V = A i? Q , 
valid for R > R m i n , where R m in is an estimate of the 
transition radius between bulge-dominated and disc- 
dominated dynamics. The analysis begins by showing 
how this model provides an extremely good descrip- 
tion of the generic behaviour of rotation curves in the 
optical disc and, furthermore, how it imposes very 
detailed correlations between the free parameters, A 
and a, of the model. 

These correlations are investigated, and shown to 



imply, via first and second-order models, a third- 
order model according to which the rotation velocity, 
V, at any radial displacement in the optical disc of 
any given spiral galaxy is given by V/Vq = (R/Ro) a , 
where Rq > R m in, and Vb are given as approximate 
functions of the galaxy's absolute magnitude and sur- 
face brightness whilst a is an unidentified function of 
other galaxy parameters - of which the most signifi- 
cant ones will be the relative proportions of the disc, 
bulge and halo mass-components. It is this latter 
function which provides the opportunity for a dark- 
matter modelling process which is independent of any 
particular dynamical theory. 

Furthermore, it is shown that the conclusion of PS 
1986, that optical-disc dynamics contain no signa- 
ture of the transition from disc-dominated dynamics 
to halo-dominated dynamics, is extremely strongly 
supported by this analysis. 



1 Introduction 

It has been know for several years - for example, Ru- 
bin, Burstein, Ford & Thonnard, hereafter RBFT - 
that optical rotation curve shapes are strongly de- 
termined by luminosity in the sense that, within 
any given morphological class and away from the in- 



ncrmost part of the rotation curve, high luminosity 
galaxies have almost flat rotation curves, whilst low 
luminosity galaxies have more 'rounded and rising' 
rotation curves. 

However, until recently all studies of the correla- 
tions existing between rotation curve kinematics and 
global galaxy properties have necessarily been con- 
fined to relatively small samples of rotation curves, 
and have therefore been subject to the corresponding 
uncertainties. This problem has been rectified by the 
publication of 900 good quality rotation curves by PS 
1995. A first analysis of this data was given by Persic 
et al (1996), hereafter PSS, in which they concluded 
that there exists a 'universal rotation curve' for spiral 
galaxies which is determined primarily by the total 
luminosity of the galaxy concerned, and the essence of 
this conclusion is confirmed here. The present analy- 
sis of the same data differs crucially from that of PSS 
in that an early decision was taken to analyse the 
data within the context of the hypothesis that rota- 
tion curves in the optical disc (that is, away from the 
dynamical effects of the bulge) could be reasonably 
described by a generalized power law, V ro t — AR a , 
where R > R m in and R m in represents an estimate of 
the transition radius from bulge-dominated to disc- 
dominated dynamics. This decision, which was mo- 
tivated by arguments which are independent of the 
matters at hand, could easily have resulted in rel- 
ative failure in the sense of providing no particular 
insights into rotation curve structure. In fact, the op- 
posite happened and it was found that the hypothesis 
imposes extremely strong correlations between rota- 
tion curve kinematics, total luminosities and surface 
brightnesses, thereby giving a detailed confirmation 
to the general thrust of the PS studies, and a strong 
confirmation to the particular idea of the universal 
rotation curve. 

In detail, the primary conclusion of this analysis 
is that the rotation curve in the optical disc of any 
given spiral galaxy with absolute magnitude M and 
surface brightness S behaves according to 
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logV 
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Ro 
-0.584 - 0.133 M - 0.000243 S, 

-3.291 - 0.208 M - 0.00292 S, 



(1) 



a rj g(d m , b m , h m ), 

where g(d mi b mi h m ) is some undetermined function 
for which the most significant parameters are prob- 
ably the disc-mass (dm), bulge-mass (b m ) and halo- 
mass (h m ) respectively. Since this model is shown 
to account for over 90% of the variation in the piv- 
otal diagram (equivalent to a regression correlation 
of 0.95), it can be considered as, at the very least, 
an extremely good approximation to the statistical 
reality, as judged for a large number of spiral discs. 

The best model that can be given for the unde- 
termined exponent, a, in terms of the data provided 
by PS 1995 is that a = 2.56 + 0.105 M, which ac- 
counts for over 29% of the variation in the (a, M) 
plot (equivalent to a regression correlation of 0.55). 
We see immediately from this model that high lumi- 
nosity galaxies are the ones which possess the flattest 
rotation curves, whilst the low luminosity galaxies are 
the ones which possess the 'rounded and rising' rota- 
tion curves - thereby confirming the RBFT result for 
a very large sample. 

A further, and very significant, consequence of the 
success of (pi) is its implied confirmation of the PS 
1986 result that optical discs appear to contain no sig- 
nature of the transition from disc-dominated dynam- 
ics to halo-dominated dynamics. PS reached this con- 
clusion after an analysis of 42 rotation curves showed 
the near-constancy, over the whole of each optical ro- 
tation curve, of the parameter V = (CI — K/2)R, for 
angular velocity CI and epicyclic frequency K . The 
connection between this result of PS and the present 
perspective is discussed in Appendix [AL 

Further evidence supporting the apparent absence 
of any kinematical transition region in the optical disc 
arises through the discovery of very strong correlation 
between R opt and R m in, where R m in is a theory- 
independent estimate of the transition region be- 
tween bulge-dominated and disc-dominated dynam- 
ics (see Appendix |b|). The existence of such a cor- 
relation implies a causal connection between R m in 
and R op t which, again, argues against the existence 
of sharply delineated transition regions between disc- 
dominated and halo-dominated dynamics. 

Finally, and of potentially greatest significance, it 
is to be noted that halo-mass only enters the sys- 



tern through the undetermined function a; this fact 
provides the possibility of a dark-matter modelling 
process which is independent of any particular dy- 
namical theory. Specifically, since, for any given 
galaxy, reasonable estimates for d m and b m can be 
obtained, then it becomes possible to model a in 
terms of these two quantities. Analysis of any system- 
atic differences arising between a-values calculated 
from the rotation curves, and the modelled a-values, 
will then, in principle, allow systematic estimates of 
any given galaxy's dark-matter content to be given 
in terms of the galaxy's disc-mass, bulge-mass and 
surface brightness. 



2 The Data 



The data given by PS is obtained from the raw Ha 
data of MFB by deprojection, folding and cosmolog- 
ical redshift correction. For any given galaxy, the 
data is presented in the form of estimated rotational 
velocities plotted against angular displacement from 
the galaxy's centre; estimated linear scales are not 
given and no data-smoothing is performed. 

The analysis proposed here requires the linear 
scales of the galaxies in the sample to be defined 
which, in turn, requires distance estimates of the sam- 
ple galaxies from our own locality. This information 
is given in the original MFB paper in the form of three 
variations of the Tully-Fisher, hereafter TF, distance 
estimate, say TF1, TF2 and TF3. TF1 is derived 
by using the TF relation for Fornax obtained from 
magnitude data; TF2 is the Malmquist-bias corrected 
form of TF1 giving TF2=1.08xTFl for all galaxies; 
TF3 is obtained by using the TF relation for For- 
nax obtained from rotation curve data. The current 
analysis is based upon TF3. 

Table 3 gives the morphological type-distribution 
in the PS 1995 data base, and shows that the great 
majority of the selected galaxies are of types 3,4,5 
and 6, with only two examples of types 0,1,2 and a 
tail of 31 examples of types 7,8,9. 



Table 3 


Galaxy 


Sample 


Type 


Size 


0,1,2 


2 


3 


306 


4,5 


177 


6 


384 


7,8,9 


31 



For the purposes of the following analysis, no distinc- 
tion is made between type classes, and so the whole 
sample of 900 rotation curves is used. However, the 
results of this analysis remain unchanged, except in 
the numerical details, when it is repeated on any of 
the subsets consisting of types {3}, {4,5} or {6}. 

The details of the data-reduction used in these 
analyses are given in Appendix H and, briefly, they 
amount to a prior- decided means of minimizing the 
effect of the bulge on the rotation-curve calculations. 

3 A Necessary Condition For 
The Power-Law Hypothesis 

Suppose the power-law hypothesis of ([j]) to be true; 
in this case, and excepting for the inevitable high 
level of noise, a plot of any given rotation curve in 
the (logi?, logF) plane would lie on a straight line. 
Since the rotation curves are not identical then, for 
the 900 rotation curves, we would obtain 900 distinct 
straight-line plots in the (logi?, logF) plane. In the 
following, we show that, in this case, the mean plot 
of the 900 separate plots must also be a straight line: 
For a given log R, the corresponding log V value for 
each of the 900 rotation curves would be given by 

log^.=a fe + 6 fe logi?, fc=1..900. 

Summing over k, and forming the mean response, we 
find 

Mean (log V ) = A + B log i?, 

so that, as stated, the mean plot is also a straight 
line. It follows that a necessary condition for opti- 
cal rotation curve data to be described by a power 
law is that the mean plot of all the rotation curves 



in the (logi?, logV) plane must be a (statistically) 
straight line. So, the strategy of the initial analysis 
is described as follows: 

• Reduce all the data in the combined sample of 
900 rotation curves to a uniform linear scale 
based on TF3; 

• Superimpose all the data of the combined sample 
into a single data set; 

• Divide the data set into bins subject to some 
convenient criteria. For present purposes, the 
bin- width was chosen as 0.057kpc subject to the 
constraint that each bin contained at least 200 
data points. In practice, the average per bin is 
about 800 data points. 

• Form the average of the data in each bin and 
plot it. 

The rationale is simply that, if the power law hypoth- 
esis is correct, then the internal noise on the data for 
the 900 separate rotation curves should largely cancel 
out over the averaging process. The result is plotted 
in Figure [l| - and it is to be emphasized that this 
plot does not represent a physical rotation curve, but 
is merely a convenient graphic tool designed to illus- 
trate a statistical point. The only discernible sense 
of any possible systematic deviation from a purely 
linear behaviour occurs in the innermost part of the 
plot; apart from this, the figure provides solid sup- 
port for the idea that a general power law is, at the 
very least, a good working approximation to H a rota- 
tion curve data, and provides the necessary rationale 
for the basic analysis of this paper. 

4 The Minimization Of Bulge- 
Dominated Dynamics In Ro- 
tation Curve Data 

Before beginning the analysis proper, we note that 
the apparent deviation from linear behaviour on the 
innermost part of the rotation curve, noted in the 
previous section, is easily understood in terms of the 
transition from bulge-dominated to disc-dominated 



dynamics. In this section, we describe the effect of 
implementing the data-reduction process described in 
Appendix H which is designed to minimize the effects 
of bulge-dominated dynamics in the present consid- 
erations . In effect, this process is a numerical pro- 
cess which removes the innermost data points on any 
given rotation curve according to how statistically 
'unusual' those points are in relation to the rest of 
the data, and its effects are described below in Fig- 
ures H & 3 and in Tables 1 & 2. 

Figure 2 plots the measured Rmin against the mea- 
sured R op t before any data reduction is performed, 
whilst Figure [| plots the calculated R m in arising 
from the data reduction process against the measured 
R op t- It is clear from these diagrams that the data 
reduction process reveals a fairly strong correlation 
between R m in and R op t- The change between the 
two diagrams is quantified in Tables 1 & 2. 



Table 1 


Ropt = bo + b\ R m in 

Before data reduction 


Predictor 


Cocff 


Std Dcv 


t-ratio 


P 


Const. 

-LMnin 


7.25 
1.62 


0.13 
0.11 


57 
14 


0.00 
0.00 


R l = 19.2% 


Table 2 


After data reduction 


Predictor 


Cocff 


Std Dev 


t-ratio 


p 


Const. 

**>min 


5.73 
1.64 


0.14 
0.07 


41 
24 


0.00 
0.00 


R z ---- 38.6% 



It is surprising to note that Table 1 clearly shows that 
Rmin in the unreduced data is already a very signifi- 
cant (if very noisy) predictor for R op t- However, after 
the data reduction process, both the significance and 
the reliability of R m in as a predictor for R op t increase 
considerably. 

If we consider the post data-reduction values of 
Rmin as strong predictors of bulge-radius (which is 
reasonable), then the strong post data- reduction cor- 
relation between R m in and R op t can be seen as a 
correlation between bulge radius and optical disc ra- 
dius. The fact that a significant correlation between 
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Figure 1: Mean Of 900 Individual Rotation Curves 
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Figure 2: Before Data Reduction 
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Figure 3: After Data Reduction 



Rmin and R op t already existed prior to the data re- 
duction can then be interpreted to indicate that the 
pre data-reduction R m in values were already signifi- 
cant (if very noisy) predictors of bulge-radius. 

There is a further significant implication of the 
strong post data-reduction correlation between R m in 
and R pt- the correlation implies a strong causal con- 
nection between the two values which, in turn, im- 
plies that any effect the halo has on R op t is necessar- 
ily balanced by a corresponding effect of the halo on 
R m in- In particular, the correlation is inconsistent 
with the existence of any strong transition region be- 
tween disc-dominated and halo-dominated dynamics 
in the optical disc - a conclusion already reached by 
PS 1986, and others and clear from Figure p]. 



5 A Basic Correlation 
Imposed By Power-Law 
Rotation Curves 

In the following, it is shown that an extremely strong 
correlation exists between the power-law parameters 



over the whole PS sample, and the correlation is de- 
tected in the following way: 

• The basic assumption is that rotation velocities 
behave as V = AR a so that log V = log A + 
a logi?; 

• For each of the 900 rotation curves in the sam- 
ple, regress log V on log R to obtain estimates of 
log A and a; 

• Plot the 900 pairs (a , log A) on a single diagram. 

The results of this exercise are shown in Figure 0, and 
shows that there exists an extremely strong negative 
(a, log A) correlation. The only obvious deviations 
from what could be called an almost perfect linear 
correlation are, firstly, a slight hint of curvature along 
the lower boundary of the distribution and, secondly, 
the fan-like behaviour of the diagram going from a 
broad spread of points at the bottom right-hand of 
the figure to a narrow neck at the top left-hand of 
the figure. 

In the following sections, we show that the struc- 
ture of Figure orders the galaxies in the sample 
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Figure 4: (a, log A) plotted for each of 900 rotation curves 



according to their absolute magnitude and surface meet at the fixed point (xo,yo) in the plane. Conse- 
brightness properties. quently, in the present context, the first-order linear 

approximation (0) implies that all rotation curves, 



6 The First-Order 70% Model 

6.1 The Basic Analysis 

The first-order model simply assumes that the 
(a, log A) relationship underlying Figure is linear, 
so that 



log A = a + b a, 



(2) 



where ao and 60 are constants, and linear regression 
shows that this model can account for 70% of the 
variation (equivalent to a regression correlation coef- 
ficient of 0.84) in Figure £|. 

6.2 Geometric Implications 

It is easily shown that when the parameters, (m,c), 
of a set of lines y = mx + c in the (x, y)-plane are 
constrained to satisfy c = j/o ~~ xqiti, for constants 
(xo,yo), then the lines themselves are constrained to 



log V = log A + a log R, 



(3) 



for which log A and a are related by (g) , intersect at 
the fixed point (— feoj a o) in the (log R, log V) plane. 
If this fixed point is denoted as (logi?o,log Vb), then 
(0) is more transparently written as 



log A = log Vo — a log Rq 
so that (||) implies 
R 



V = V 



(4) 



(5) 



Assuming the model (0), a linear regression on the 
data of Figure || subsequently gives (log i? , log V ) w 
(0.84,2.22). Whilst the idea of a single intersec- 
tion point for all rotation curves in the sample seems 
rather extreme, it is unambiguously deduced as a con- 
sequence of the first-order modelling of Figure - and 
is therefore to be seen as a first-order approximation 
of the reality. 



In the following sections, we give two powerful illus- 
trations of the reality of the idea that optical rotation 
curves can be considered, in the first approximation, 
to converge in a very small region of the (log R, log V) 
plane. 

6.3 The Intersection Frequency 
Diagram 

A direct illustration of rotation curve convergence is 
given as follows: we have calculated the coordinates 
of intersection between all possible pairs of 200 ran- 
domly sampled rotation curves drawn from the full 
sample, and have plotted the frequencies of intersec- 
tion along the logi? and logV axes respectively in 
Figures |] & [| It is clear from these figures that there 
is a very sharp peak of intersection points in both di- 
agrams. Using the mean points in each diagram, to- 
gether, they are consistent with a peak intersection 
point at (log R, log V) sw (1.02, 2.24). 

6.4 A Geometric Illustration Of Rota- 
tion Curve Convergence 

A second, and dramatic, illustration of rotation curve 
convergence can be formed from the realization that, 
if the rotation curves really do converge on a single 
fixed point in the (log R, log V) plane, then all of the 
rotation curves will transform into each other under 
rotations about the fixed point, (log.R ,log Vb), in 
this plane. 

This idea can be tested in the following way: sup- 
pose that an arbitrarily chosen straight line passing 
through the point (log Rq, log Vb) is defined as a stan- 
dard 'reference line'. Since, according to the first- 
order model interpretation of Figure [|, all the ro- 
tation curves in the sample pass through the point 
(logi?o, log Vb), then every rotation curve in the sam- 
ple can be transformed into the standard reference 
line by a simple bulk-rotation about (log.Ro: log Vb). 
Since, according to this idea, the rotation curves are 
reduced to equivalence by the rotation, then the pro- 
cess of forming an 'average rotation curve' from the 
set of rotated such curves should greatly reduce the 
internal noise associated with the individual rotation 
curves, and we would expect the resulting average 



curve to be a very close fit to the standard reference 
line, referred to above. 

In the following analysis, which is performed for all 
900 galaxies in the sample, the technique described 
in Appendix O is used to refine the histogram esti- 
mate of (1.02,2.24), given in the last section, to the 
optimal value for the fixed point as (log.Ro: log Vb) = 
(1.06, 2.29). The solid line in Figure ^is the standard 
reference line which is defined to pass through the 
fixed point (log.Ro, logVb) = (1.06,2.29), and fixed 
arbitrarily so that log A — 1.5. The results are not 
sensitive to the choice of this latter parameter. The 
circles in the figure show the results of rotating the 
900 rotation curves of the sample about the fixed 
point so that they coincide in a least-square sense 
with the standard reference line, and then averaging 
over all the rotation curves. Each circle represents 
an average of approximately 800 separate measure- 
ments. 

When Figure [7] is compared with Figure [l], we see 
that virtually all of the scatter present within the lat- 
ter figure is eliminated in Figure PL thereby provid- 
ing the strongest possible evidence for the idea of the 
equivalence of rotation curves with respect to rota- 
tion about particular fixed points in the (log R, log V) 
plane. 



6.5 A Comparison of Three 

(log.Ro, logVo) Estimates 



As we pointed out in §6.1, the first-order approxima- 
tion of the (a, log A) plot of Figure [| implies directly 
that all H a rotation curves pass through a single fixed 
point in the (log R, log V) plane, and we have three 
distinct ways of estimating the position of this fixed 
point. The first way was by a direct linear regres- 
sion on the data of Figure 0, the second way was by 
estimating the peaks of the intersection histograms 
of Figures ra and ra respectively, whilst the third way 
was by the minimization procedure described in § 3.4 
the results of which are shown in Figure pj. The three 
estimates are shown together in Table 4. 
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Figure 5: Frequency of Intersections on logi? axis 
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Figure 6: Frequency of Intersections on logV^ axis 
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Figure 7: Types 0..9 spirals; Averaged rotated data 



Table 4 


Method 


logi? 


logVo 


Regression 
Mean Values 
Minimization 


0.84 
1.02 
1.06 


2.22 
2.24 
2.29 



correlation coefficient of 0.92) in the pivotal diagram, 
Figure ||. 



There is clearly a high degree of consistency between 
the three methods which serves to emphasize the re- 
ality of the first-order effect. 



7 The Second-Order 85% 
Model 

7.1 The Basic Analysis 

Whilst, so far, Figure Hhas been analysed on the ba- 
sis of a first-order linear model, there is the possibility 
that systematic internal structure exists in Figure 0. 
In this section, we show how the existence of this 
internal structure is strongly supported on the data, 
and leads to a second-order model which accounts for 
over 85% of the variation (equivalent to a regression 



Table 5 


(a, log A) data averaged in 


magnitude-limited quartiles 


Sample 


Mean 


Mean 


Mean 


Size 


abs(mag) 


a 


log A 


224 


-22.20 


0.22 


2.15 


227 


-21.16 


0.34 


1.96 


224 


-20.27 


0.47 


1.79 


225 


-18.80 


0.58 


1.65 



We begin by taking the 900 pairs of (a, log A) data 
which form the basis of Figure ^ and partition them 
as near as possible into quartiles according to ab- 
solute magnitude; thus the data is partitioned into 
the brightest 25% of the sample, the next brightest 
25% and so on. For each quartile, we then form the 
mean values of the absolute magnitude, log A and a 
respectively, and the results of this exercise are listed 
in Table 5 and plotted in Figure |l The solid line in 
the figure is the linear regression line fitted to the four 
data points, and the whole diagram can be considered 
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Figure 8: (a, log A) data averaged in magnitude-limited quartiles 



as a condensed representation of the data plotted in 
Figure ^ showing the variation of absolute magnitude 
through that data. It is clear from Table 5 and Fig- 
ure |8J that a and log A are each very strongly corre- 
lated with the absolute magnitude, thereby justifying 
a more detailed investigation. 



7.2 The Refined Analysis 

We begin by plotting the (a, log A) data for each 
of these magnitude-limited quartiles and displaying 
these plots in Figure @. Viewed collectively, the four 
diagrams of this figure display three significant fea- 
tures: 

• Firstly, the diagram corresponding to the 
brightest magnitude-limited quartile shows that 
(a, log A) data exhibits almost scatter-free lin- 
ear behaviour, whilst the data of the remaining 
quartiles exhibits increasing scatter as absolute 
magnitude increases. 

• Secondly, apart from this scatter, there appears 
to be a systematic variation of the data through 



the four diagrams when they are viewed in order 
of increasing absolute magnitude. 

• Thirdly, the scatter evident in Figure H is seen to 
be a function of increasing absolute magnitude 
which, at face value, would suggest it is a func- 
tion of increasing measurement uncertainties at 
the dimmer end of the sample. However, there 
is a feature which suggests that the situation is 
not this straightforward: specifically, the lower 
boundary of points in each diagram of Figure 
M is very sharply delineated whilst it is the up- 
per boundary which becomes increasingly diffuse 
in the higher magnitude quartiles. If measure- 
ment uncertainty is the basic cause of scatter in 
these diagrams, then it is not immediately obvi- 
ous why the lower boundary of points in every 
quartile is sharply delineated. 

The reality of a systematic variation with increasing 
magnitude in the diagrams of Figure fl is confirmed 
by forming the linear regression of log A on a for the 
data of each of these diagrams so that, for each dia- 
gram, there is a linear model analogous to (0). The 
results are listed in Table 6 and plotted in Figure HOl 
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Figure 9: (a, log A) plotted for each magnitude-limited quartile 



The figure shows very clearly that the four regression 
lines meet in a very small neighbourhood of the point 
(0.9, 1.5) of the (a, log A) plane. 



Table 6 


log A = log Vo — a log i?o 


Sample 
Size 


logVb 


log R 


Mean M 
of quartile 


224 
227 
224 
225 


2.35 
2.23 
2.10 
1.89 


0.90 
0.80 
0.65 

0.41 


-22.20 
-21.16 
-20.27 
-18.80 



We have already noted that, when a set of lines 
y = mx + c in the (x, y)-plane are constrained to 
meet at a fixed point, (xo,yo) say, then the param- 
eters (m, c) of the lines are constrained to satisfy 
c = yQ — XQjn. In the present case, representing any of 
the lines in Figure |ltj| as log A — log V — a logR , this 
means that the parameters (logi?0j log^o) & r e con- 
strained to satisfy logVo ~ 1.5 + 0.9 logi?o- Since, in 
Figure nfl, the only thing which varies between the 
lines is the absolute magnitude, then we can expect 



logVb and logi?o to be individual functions of ab- 
solute magnitude and this is confirmed by Table 6. 
Linear regression of the data of this table gives 



log V = 
logi? = 



-0.67- 0.14 M, 
-2.34 -0.15 M 



(6) 



from which we get log Vb = 1.51 + 0.93 log.Ro, which 
confirms the result obtained directly from Figure [tOl 
We now note that the constants and coefficients in 
rtq) can be 'fine tuned' in the following way: If we 
substitute equations (pi) into (H) we obtain the rela- 
tionship 

\ogA = logVb -a logi? ->■ 
log A = -0.67- 0.14 M 

+ 2.34 a + 0.15 aM, (7) 

which shows directly that the model (||) is equiva- 
lent to a model for log A expressed in terms of M, a 
and a M. Given this model for log A, its coefficients 
can be optimised in a least-square sense by regressing 
log A directly on the predictors M, a and aM over 
the whole data set. This process gives 

log A = -0.56 -0.13 M 
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Line A: M = -22.20 
LineB: M = -21.16 
Line C: M = -20.27 
LineD: M = -18.80 
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alpha 

Figure 10: Regressions on (a, log A) data in magnitude-limited quartiles 



log A 
logVb 

logi? 



2.27a + 0.14 a M, -> 
log Vq — a log R , where 
-0.56- 0.13 M, 
-2.27- 0.14 M. 



with very strong statistics on the constant and coef- 
ficients, and shows that this regression accounts for 
over 85% of the variation (equivalent to a regression 
correlation coefficient of 0.92) in the pivotal diagram, 
which is Figure 0. A comparison with (pi) shows only 
minor adjustments to the constants and coefficients 
of the (log.Ro, fogVo) models. Consequently, and us- 
ing (||), the optimised second-order model gives the 
universal rotation curve as 



the second-order model (g) which depends only on 
absolute magnitudes. However, the PS data-base for 
the 900 rotation curves also provides estimates to the 
optical radius, R op t, which, with the absolute mag- 
nitudes, allows us to estimate the surface brightness 
of each galaxy in the sample; since this represents in- 
dependent information, it seems sensible to consider 
the effect of its inclusion in the model and, for present 
purposes, we define it to be in solar luminosities per 
square parsec so that 

'2.5 5 - A/N 



S 



irR 2 



V_ f R_ 

log^o = -0.56- 0.13 M, 
logR = -2.27 -0. 14 M. 



(8) 



8 The Third-Order 90% Model 

It has been shown that over 85% of the variation in 
the pivotal diagram, Figure 0, is accounted for by 



The second-order model, (]8j), was arrived at by 
initially noting that Table 6 implied the relations 
(|6|) which, in conjunction with the basic relationship 
log A = logVb — a logi?o, allowed the deduction of 
the structure 

log A = ao + ai a + ai M + 03 a M 

for the log A model. If (||) is to be modified by the 
introduction of S, so that 



log^o 
log-Ro = 



-b 
cq- 



-b x M + b 2 S 
c\ M + c 2 S, 
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then the log A model will have the structure 



8.1 A Secondary Test Of The Model 



log A = b + biM + b 2 S 

— cq a — c\ a M — C2 a S. 



(9) 



The results of regressing log A on the predictors M, 
S, a, a M and a S are given in Table 7. 



Table 7 


log A = 


bo + b\ M + 62 S — Cq a — c\ 


en M — ci a S 


Prcd 


Coeff 

xlO 4 


Std Dev 
xlO 4 


t-ratio 


P 


Const. 


-5840 


890 


-7 


0.00 


M 


-1330 


44 


-30 


0.00 


S 


-2.434 


0.77 


-3 


0.00 


a 


32910 


1600 


21 


0.00 


aM 


2080 


84 


25 


0.00 


a S 


29.23 


2 


14 


0.00 


R 2 = 90.3% 



This model now accounts for marginally over 90% of 
the variation (equivalent to a regression correlation 
coefficient of 0.95) in the pivotal diagram, Figure |], 
and we see from the t-ratios that the surface bright- 
ness is most certainly a significant component of the 
model. Interpreting the table according to the model 

log A = log Vo - a log R , 
log V = b + h M + b 2 S, 
log.Ro = co + ciM + c 2 S, 



we arrive, finally, at the third-order model 
V 

logV = -0.584- 0.133 M- 0.000243 S, 
\ogRo = -3.291- 0.208 M- 0.00292 S. 



-Y 

RoJ 



(10) 



Referring to Table 7 again, we see that S is a very 
strong predictor for logi?o> whilst its effect on log Vo, 
whilst significant, is much less so. This conclusion is 
supported by numerical experimentation which shows 
that it can be omitted from log Vo with hardly any 
effect on the model. However, Table 7 shows that it 
probably is a real predictor for logVb, and so it is 
retained in the model. 



A secondary test of the validity of the third-order 
model can be given as follows: For each galaxy 
we compute (Rq,Vq) according to ([!(]), form the 
scaled profile (R/R , V/Vq) and regress log(V"/Vb) on 
log(i?/J?o)- If ( |l0l ) is good enough, then the regres- 
sion constant should be statistically zero. The dis- 
tribution of the 900 regression constants is shown in 
Figure O, for which a one-parameter t-test gives the 
95% confidence interval for the mean of the distri- 
bution as (-0.003, +0.007) whilst the 99.999% con- 
fidence interval is given as (—0.010,-1-0.014). Thus, 
for all practical purposes, the regression constant can 
be considered as statistically zero, thereby confirming 
the quality of the model. 

8.2 The Brightest Quartile 

Finally, as a means of demonstrating that the (al- 
ready small) scatter about the zero point evident in 
Figure |llj is a function of uncertainty in the data, 
rather than due to a poor fit of the basic power-law 
model, we apply the third-order model directly to 
the 224 rotation curves of the brightest quartile, for 
which the (a, log ^4) data is plotted in the first di- 
agram of Figure |9[ For this data, the distribution 
of the 224 regression constants arising from the lin- 
ear regression of log(V A /Vb) data on \og(R/Ro) data 
is given in Figure |l2|. A direct visual comparison 
with Figure [ll] shows immediately that the scatter of 
points about the zero point is very much reduced in 
the brightest quartile data, thereby supporting the 
hypothesis that the scatter which exists in Figures 
nil and O is due purely to uncertainty in the data, 
and not at all due to any lack of fit of the general 
power-law model. 

8.3 The Surface-Brightness 
Cut-Off 

It has already been noted that the lower bound- 
ary of points in each of the magnitude-limited quar- 
tile plots of (a, log A) shown in Figure S is very 
sharply delineated, whilst the corresponding upper 
boundaries become increasingly diffuse as magnitude 
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Mean Log(A) = 0.002 
Std Dev Log(A) = 0.081 
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Figure 11: Distribution of 900 regression constants 
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Figure 12: Brightest quartile regression constants 



15 



-23.55<M<-21.58 



-21.59 <M< -20.72 



2.5 



2 



1.5 





S 


=300 

M=-23.54 


s=o V 




o \. 

°o X 


M=-21.59 







0.5 1 1.5 

-20.73<M<-19.74 



2.5 






S=190 


2 




^Oq 


D 

jj.8 M=-20.72 


1.5 


S=0 


o 






M= 


19.75 


°\ o 



0.5 1 

alpha 



1.5 



2.5 



2 



1.5 







S=250 






>q M=-21.58 


S=0 




a, oo\ 


M=- 


20.73 


\ 



0.5 1 1.5 

-19.74 <M< -16.07 




0.5 1 

alpha 



Figure 13: Surface- brightness boundaries shown for each magnitude- limited quartile 



increases. We shall show that this phenomenon 
arises because the lower boundary of points in each 
of this diagrams effectively defines a zero surface- 
brightness boundary and therefore represents a real 
physical cut-off boundary. By contrast, the upper 
boundaries are closely associated with maximum sur- 
face brightnesses in each of the quartiles, and so 
the magnitude-dependent scatter observable through 
the four magnitude- limited quartiles almost certainly 
arises from observational uncertainties. This is shown 
in the following manner: for each of the quartiles, 
which is magnitude-limited in its own specific range 
{M m i n < M < M max ), we determine the surface- 
brightness limits, (0 < S < Smax), which contain 
95% of the sample (to eliminate the effects of long 
tails). Then, using the third-order model, log A = 
F(a, M, S), we find that the distribution of points 
in each of the quartiles is closely bounded by the lines 



log ,4 = F(a, M, 
log A = F(a, M, 



,0), 



mm: ^max 



). 



in the (a, log A) plane for each of the magnitude- 
limited quartiles. The results are shown in Fig- 



ure 03, and the plots in the three brightest quar- 
tiles show very clearly that the lower boundaries of 
points in these quartiles coincide very closely with the 
(S = 0, M = M max ) boundary in the plane, whilst 
the diffuse upper boundaries of points coincide very 



closely with the (S 
in the plane. 



Smax, M = M mm ) boundary 



The only partial exception to this rule is the plot 
in the dimmest magnitude-limited quartile. Here, the 
upper (very diffuse) boundary of points still coincides 
with the {S — S max , M — M m i n ) boundary; how- 
ever, it was found that the (S = 0, M = M max ) 
boundary is not even close to the lower boundary of 
points in the plot but, instead, we find that the lower- 
boundary of points follows closely the (S = 0, M = 
Mmin) boundary The important point is that the 
lower boundary of points is still a S — boundary. 
The obvious explanation for this discrepancy is that 
objects of very low surface brightness and very high 
magnitude are observationally excluded by selection 
effects, and therefore do not exist in the dimmest 
quartile sample whilst, simultaneously, the same se- 
lection effects will also ensure that the low surface- 
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brightness objects in the dimmest quartile will tend 
to be the low magnitude objects. 

To summarize, the sharply delineated lower bound- 
aries in the magnitude-limited quartiles of Figures 
or Ol define a physical 5 = cut-off boundary, whilst 
the increasingly diffuse upper boundaries probably 
arise as a function of the increasing observational un- 
certainties associated with increasing magnitudes. 

9 A Comparative Test Of 
The First, Second And 
Third-Order Models 



The three models share the basic structure 
V f R 



(11) 



and are distinguished only by the models 
assumed for (i?o, Vo). These three mod- 

els are defined, respectively, according to: 
The First-Order Model 

logi? = 1-02, (12) 

logVb = 2.29; 

The Second-Order Model 

log.Ro = -2.27- 0.14 M, (13) 

logVb = -0.56- 0.13 M; 

The Third-Order Model 

logi?o = -3.291 - 0.208 M- 0.00292 5,(14) 
logVo = -0.584- 0.133 M- 0.000243 S. 

As a final comparative test of these three models 
we note that, assuming each to be 'correct', then 
a linear regression of log(V/Vo) data on log(R/Ro) 
data should give rise to a statistically zero constant 
terms, indicating that the regression lines (tend to) 
pass through the origin. The distribution of these 
constant terms for each of the 900 rotation curves is 
given, for each of the three models, in Figure H4. To 
make the visual comparison easy, the three diagrams 
in the figure are each drawn to the same scale, and 
it is clear that the process of successively refining 



the model has the effect of tightening the distribu- 
tion of log A about the zero-point. Specifically, the 
first-order model is comparatively very poor, whilst 
the third-order model is visually significantly better 
than the second-order model. 



Table 8 


One-parameter t-test for confidence 
limits of Mean(log A) 


Model 


95% Interval 


99.999% Interval 


1st- Order 
2nd-Order 
3rd-Order 


0.009,0.029 
-0.006,0.007 
-0.003,0.007 


-0.003, 0.041 
-0.014,0.015 
-0.010,0.014 



The results of a one-parameter t-test for limits on 
the position of the mean value of log^4 for each of 
the models are given in Table 8. Apart from showing 
that the log A distribution for the lst-order model is 
slightly asymmetric with respect to the zero point, 
the table shows that the 2nd and 3rd-order mod- 
els place very tight 95% and 99.999% confidence in- 
tervals on the mean of log A about the zero point. 
Thus, for all practical purposes, the regression con- 
stant computed over the 900 rotation curves is statis- 
tically zero in both the second and third-order mod- 
els, so that these models receive the strongest possible 
support from the data. 

10 Equivalence Classes Of Ro- 
tation Curves 

As already discussed in detail in §KJ, the first-order 
model, equation ( |ll| ) with (|1^), effectively defines the 
set of all rotation curves in the sample as a single class 
within which all the rotation curves pass through 
the single point (logi?o, logVb) in the (logi?, logV) 
plane. 

By contrast, the second-order model, equation (|l|) 
with (^3[), replaces the single point, (logi?o, logVb), 
of the first-order model by a class of such points 
parametrized by M. All rotation curves in the class 
specified by M then pass through the single point, 
(log Rq, log Vo), defined by @. 

The third-order model, equation ( pi] ) with (O), is 
a refinement of the second-order model in which the 
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Figure 14: Comparative test of the three models 



class to which any given rotation curve belongs is de- 
fined by specifying M and S, the surface brightness. 
All rotation curves in the class specified by (M, S) 
then pass through the single point, (log-Rch logVo), 
defined by (£[§). 

If we make the (reasonable) assumption that all 
of the systematic variation in the pivotal diagram, 
Figure 0, has been accounted for by the 90% third- 
order model then we can realistically assume that, 
although the functional form of the model is proba- 
bly an approximation to some ideal, there are no fur- 
ther physical parameters involved in the specification 
of (logi?o, logVb). It will then follow that the rota- 
tion curves of all galaxies of having the same absolute 
magnitude and surfaces brightness will pass through 
a single point, (log.Ro, logVo), in the (logi?, logV) 
plane. Consequently, all such rotation curves will be 
equivalent to within a rotation through the point, and 
can therefore be considered to define an equivalence 
class of rotation curves. 



11 The Universal Rotation 
Curve 

Using the third-order model, it is shown that the uni- 
versal rotation curve for spiral galaxies can reason- 
ably be assumed to have the general structure 



(15) 



V 
V 


UoJ 


( 


log-Ro = 


= -3.291 -0.208 M- 


-0.00292 5", 


log Vo = 


= -0.584- 0.133 M- 


- 0.000243 S, 


a = 


CJ\&rni Urn-, ft"m) 





where g(d m , b mi h m ) is some undetermined function 
for which the most significant parameters are proba- 
bly the disc-mass, bulge-mass and halo-mass respec- 
tively of the spiral galaxy concerned. 

11.1 The General Situation 

Reference to (0) with (|l|), (|l|) and |L4|) shows that, 
once a model for a in terms of the physical proper- 
ties of spiral galaxies is obtained, then (O) with any 
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one of (|12[), ( |13[ ) or ( |14| ) can be considered to rep- 
resent a universal rotation curve. In the following, 
we consider the problem within the context of the 
third-order model. 

As we have already noted, the issue depends en- 
tirely upon the function a, and what galaxy param- 
eters it is a function of. 



the model for a described in Table 9 so that, specifi- 
cally: 



Table 9 


a = b a + bi M 


Predictor 


Coeff 


Std Dev 


t-ratio 


P 


Const. 

M 


2.560 
0.105 


0.112 
0.005 


23 
19 


0.00 
0.00 


R 2 = 29.2% 



A least-square linear modelling of a in terms of the 
available parameters, M and S, gave the best model 
as a = —0.583 + 0.105 M, with no evidence of any 
significant dependency on S. The details of the best 
model are given in Table 9 which shows the depen- 
dency of a on M to be extremely significant, but 
which also shows that the model accounts for only 
29% of the total variation in the (a, M) plot (equiv- 
alent to a regression correlation of approximately 
0.55). To test the possibility that the scatter in this 
latter plot was, perhaps, due to noisy data, we tried 
modelling on the brightest 25% and 50% of the data, 
but found no evidence at all to support this possi- 
bility. The only reasonable conclusion was that the 
unaccounted-for 71% variation in the (a, M) plot 
arises through the effects of unaccounted-for addi- 
tional physical parameters. 

The major unincluded determinants of disc dynam- 
ics are the relative amounts of disc-mass, halo-mass 
and bulge-mass. Consequently defining d m = disc- 
mass, b m = bulge-mass and h m = halo-mass, and 
noting that M is effectively a measure of the to- 
tal visible mass, then we can reasonably assert that 
a s=s g(d m , b m , h m ). 

11.2 A Simple Approximation 

Finally it is interesting to note how the foregoing 
analysis allows the construction of a very simply uni- 
versal rotation curve, based on absolute magnitude 
alone. This is given by (O) and (Of) together with 



V 

~% 
log-Ro 

logVb 

Q 



R 

Ro , 

-2.27- 0.14 M, 
-0.56- 0.13 M, 
2.560 + 0.105 M, 



(16) 



which can be considered as a formal refinement of 
the information contained in Figure p[ This latter 
approximation confirms the conclusions of PS 1996 
which are that the universal rotation curve is sub- 
stantially determined by absolute magnitude alone, 
especially for the brightest galaxies. 



12 Theory-Independent 
Matter Models 



Dark- 



PSS, used essentially the same data set as that anal- 
ysed here to show that whilst, at high luminosities, 
there is only a slight discrepancy between observed 
rotation curves and those predicted on the basis of the 
observed luminous matter distributions, the discrep- 
ancy is far more serious at low luminosities. This con- 
clusion is already strongly supported by the simple 
model given at (nq) which shows how the power-law 
exponent, a, increases strongly with absolute magni- 
tude. Since, (relatively) large values of a correspond 
directly to the most steeply rising rotation curves, 
and since it is precisely these kinds of rotation curves 
that require the largest amounts of dark matter for 
their explanation, according to the virial theorem, 
then the findings of the present analysis are in direct 
accord with those of PSS. 

In this way, PSS have already used the data to 
show that the dark-matter component of spiral galax- 
ies does not seem to have the form of an uncon- 
strained and independent structure added onto the 
visible structure of spiral galaxies, but appears to 
make its presence felt in some kind of systematic way 
which is strongly inversely correlated with the visi- 
ble component of the galaxy structure. This effect, 
together with the detailed considerations of the fore- 
going sections, allows the construction of dark- matter 
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models which are independent of any particular dy- 
namical theory and therefore, as a side effect, provide 
strong tests of such theories. 

Specifically, referring to the final form of the uni- 
versal rotation curve, given at ( |l5| ) , we see how dark- 
matter can only make its presence felt through the 
exponent, a. Given that reasonable light-based esti- 
mates can be made of disc-mass (d m ) and bulge-mass 
(b m ) for spiral galaxies, then it is (in principle) a sim- 
ple process to form empirical models of a directly in 
terms d m and b m . Once such models are formed for 
a reasonably sized data-base of spiral galaxies, then 
it becomes a straightforward process to investigate 
the structure of any systematic deviations of predic- 
tions based on the model a « g{d m , b m ) from the 
a- values calculated directly from the rotation curves. 
The presence of any such systematic deviations can 
then reasonably be considered due to the presence of 
a dark-matter halo, so that, subsequently, estimates 
for the amount of dark-matter present for any given 
galaxy can be given in terms of the directly measur- 
able properties of the galaxy concerned. 

13 Conclusions 

By demonstrating how completely the power-law 
model for rotation curves resolves rotation curve 
data, magnitude data and surface brightness data 
over a large data-base, this paper: 

• substantially refines the well-known correlations 
that exist between the rotational kinematics of 
spiral galaxies and various of their non-kinematic 
properties, such as absolute magnitude; 

• shows how the power-law rotation curve is, al- 
most certainly, exact for idealized discs (that 
is, for discs without the irregularities inevitably 
present in real optical discs). 

Additionally, at its lowest approximation, the analy- 
sis confirms the conclusion of PSS (1996) that the ro- 
tation curve of any given galaxy is largely determined 
by the absolute luminosity of the galaxy concerned, 
and supports their conclusion that absolute lumi- 
nosity is inversely correlated with dark mass. Such 



an inverse correlation obviously provides a means of 
dark-matter modelling, but the detailed nature of the 
present analysis allows a further step in providing a 
means for the generation of statistical dark-matter 
models which are independent of any particular dy- 
namical theory. Apart from the objective desirability 
of such models, they have an obvious potential in dis- 
tinguishing between competing dynamical theories. 

Furthermore, in showing that optical discs contain 
no signature of the transition from disc-dominated 
dynamics to halo-dominated dynamics, and in show- 
ing the existence of a strong correlation between 
Rmin and R op t, the work provides very strong sup- 
port for the conclusion of PS (1986) that the disc- 
mass distribution and the halo-mass distribution are 
very strongly correlated, and suggests that the bulge- 
mass distribution should also be factored into this 
correlation. 



A The Constancy of V 

PS (1986) showed that the parameter 
V 



n-£)n. 



where Q is angular velocity, and K, the epicyclic fre- 
quency, is defined by 



K 2 



1 d 
R3~dR 



(V 2 R 3 ) 



is almost constant over the larger part of any given 
optical disc for a sample of 42 spiral galaxies, and 
they deduce from this constancy that rotation curves 
contain no signature indicating transition from disc- 
dominated dynamics to halo-dominated dynamics. In 
the present case, the same conclusion can be drawn 
from the fact that the kinematics over the whole opti- 
cal disc are amenable to description by a single power- 
law prescription. Such a power-law also provides an 
extra insight into the constancy of the V parameter: 
Using the general power-law prescription V — AR a 
in the above we find, after some algebra, 



dV_ 
~d~R 



AR 



a-0.5 



R 



V2, 



20 



From this, we see that, whenever whenever a rotation 
curve is such that a < —0.5, then the correspond- 
ing dV/dR rapidly becomes small for increasing R - 
and this corresponds to a nearly constant V . Ref- 
erence to Figure @ shows that this conditions holds 
for all galaxies brighter than M rs —20 and, of the 
PS (1986) sample of 42 galaxies, 37 are brighter than 
M w —20 so that a < —0.5 generally, indicating very 
flat V functions, as actually shown by PS. 



B Data Reduction 

A basic assumption is that optical rotation curves 
can be described in terms of a simple power law, and 
the validity of this assumption as, at the very least, 
a good approximation to the reality over the approx- 
imate range (1, 20)kpc has been demonstrated in §||. 

However, spiral galaxies can be considered to con- 
sist of a spherical bulge surrounded by a flat disk, 
and we can be reasonably certain that the bulge and 
the disk represent distinct dynamical regimes. In 
practice therefore, the structure of the inner parts 
of rotation curves will be determined one dynami- 
cal regime, whilst the structure of the outer parts 
of the rotation curves will be determined by another. 
Since it is not reasonable to suppose that a single sim- 
ple power law can bridge the two dynamical regimes, 
there is a need to find a way of eliminating the ef- 
fect of the bulge dynamics from the rotation curve 
data. In order to ensure that no subjective bias is 
imposed on the data during such a data reduction 
process, the process adopted must be automatically 
applicable according to strictly predefined rules. The 
adopted procedure is described below. 

The original analysis was done using the commer- 
cial statistics package, Minitab; when regressions are 
performed in Minitab, it automatically flags obser- 
vations which have an unusual predictor or if they 
have an unusual response. In both cases 'unusual' 
is defined by internal parameters of the software. It 
was decided, in advance of any data processing, that 
when a Minitab regression flagged the innermost ob- 
servation on any given rotation curve (that is, the 
one most likely affected by the central bulge) as 'un- 
usual' then that observation would be deleted from 



the analysis, and the regression repeated. If, on the 
repeated regression, the new innermost observation 
was flagged as unusual, then it too would be deleted 
from the analysis. This process was repeated until 
the innermost observation remained unflagged after 
the regression process. 

There are 900 rotation curves with a total of 19183 
observations recorded with a non-zero radial coordi- 
nate (those with a zero radial coordinate were nec- 
essarily deleted because data was transformed into 
logarithmic form). Of these 19183 observations, the 
'innermost deletion' strategy led to the exclusion of 
2264 observations, which is 11.8% of the sample. This 
was considered to be an acceptable attrition rate. 

C Determination of Fixed 
Points, (log.Ro, logV ) 

In §4 and §5 it was determined that the form of Fig- 
ure [| implied the first-order hypothesis that all ro- 
tation curves in the sample passed through a fixed 
point, (logi?o,logVo), in the (logi?, logF) plane. 
The problem is to determine the position of the fixed 
point (log-Ro, logVb). In the following, we describe 
how this can be done using a minimization procedure. 

The hypothesis implies that all 900 rotation curves 
in the sample should be equivalent to within a ro- 
tation about (log-Ro, logVo); in this case it should 
be possible to superimpose all such rotations curves 
upon each other by means of a simple appropriate 
rotation. However, the very noisy nature of rotation 
curve data means that, at best, such a process of su- 
perimposition can only be defined in some averaged 
statistical way. 

In practice, what is done can be characterized as 
follows: 

• Define an initial guess for the fixed point 
(log.Ro, logVb); this is done by the simple ex- 
pedient of identifying the peaks of the rotation- 
curve intersection frequency distributions show 
in Figures || & f| and they are found to be 
logi?« 1.02 and log Vw 2.24. 

• Define a standard reference curve passing 
through the estimated fixed point; typically, this 
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is done by making the arbitrary choice log A = 
1.5 and then using (0) to determine a. 

• Rotate the data corresponding to each in- 
dividual rotation curve about this estimated 
(log-Ro; log Vq) until it coincides in a least-square 
sense with the standard reference curve. 

• Keep a cumulative sum of all the least-square 
residuals arising from the foregoing rotation op- 
erations performed on the individual rotation 
curves. 

• Minimize this cumulative sum of least- 
square residuals with respect to variations 
in the estimated position of the fixed point, 
(log-Ro, logFo). 

Note that the minimization processes referred to 
above are minimizing functions defined on noisy data. 
It is routinely recommended that such problems are 
solved using the Simplex method, which is extremely 
robust - although slow. The original implementation 
of the method is given by Nelder & Mead (1965). 
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